Destruction of Neel order and local spin spirals in insulating La2-xSrj:Cu04 
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Starting from the t-J model, we derive an efltective field theory describing the spin dynamics in the 
insulating phase of La2-a;Sra;Cu04, x < 0.055, at low temperature. Using Monte Carlo simulations, 
we show that the destruction of Neel order is driven by the single-hole localization length k. A phase 
transition at 2% doping is consistent with the value of k known from the variable range hopping 
conductivity. The static spin structure factor obtained in our calculations is in perfect agreement 
with neutron scattering data over the whole range of doping. We also demonstrate that topological 
defects (spin vortex-antivortex pairs) are an intrinsic property of the spin-glass ground state. 
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Introduction. One of the most intriguing properties 
of La2-2;Sr2;Cu04 (LSCO) is the static incommensurate 
magnetic ordering at low temperature observed in elas- 
tic neutron scattering experiments. This ordering man- 
ifests itself as a scattering peak shifted with respect to 
the antiferromagnetic position: Q = Qaf + <5Q, where 
Qaf = (tt, ±7r), setting the lattice spacing a — 1 and 
using the tetragonal coordinate system. The incommen- 
surate ordering is a generic feature of LSCO. According 
to experiments in the Neel phase, x < 0.02, the incom- 
mensurability is almost doping independent and directed 
along the orthorhombic b axis, 5Q « ±0.02V2(7r,-7r), 
see Ref. In the spin-glass phase (0.02 < x < 0.055), 
the shift is again directed along the orthorhombic b axis, 
but scales linearly with doping, SQ ~ ±\/2x{n,—n), see 
Refs. 0,0,11. Finally, in the undcrdoped superconduct- 
ing region (0.055 ^ x < 0.12), the shift scales linearly 
with doping and is directed along the crystal axes of the 
tetragonal lattice, 6Q ~ 2x(±7r, 0) or SQ « 2a;(0,±7r), 
see Ref.H 

These observations caused a renewal of theoretical 
interest in the idea of spin spirals in cuprates, both, 
from t he p henomenological 0, 0, and the micro- 
scopic |lOl im ll^ point of view. The jump of the spi- 
ral direction at the insulator-superconductor transition 
(x sa 0.055) was explained in Ref. and a mechanism 
for the pinning of the spiral direction to the orthorhom- 
bic b axis in the insulating phases (Neel and spin-glass 
regime) was suggested in Ref. In this letter, we formu- 
late an effective field theory that describes the spin-glass 
phase of LCSO at low temperature. We study the ground 
state, calculate the spin structure factor, and analyze the 
nature of the quantum phase transition from the Neel to 
the spin-glass state. 

Model. In insulating LSCO, holes are trapped by 
Sr ions and form hydrogen-like bound states ip (r) — 
(i") — ^•\/2/7rKe~'"', where is a two-component 
spinor describing the pseudospin and k « 0.4 is the 
inverse localization leng th We refer to these bound 
states as impurities |l4l | . The temperature corresponding 



to the binding energy is about Tf, ^ lOOK llj; hence at 
T <^Tb the charge degrees of freedom are frozen. Due to 
the small orthorhombic distortion of LSCO, impurities 
in the Sb = (7r/2,— 7r/2) hole pocket have lower energy 
than impurities in the vicinity of Sa = (7r/2, 7r/2). From 
the estimation of the pinning energy Il2l |. we conclude 
that for T <^ Tp 30K, impurities are pinned to the 
Sb valley and the spirals are directed along the b axis, 
as explained in Ref. 11. Since we consider T <^ Tp, the 
relevant degrees of freedom are the spins of the Cu ions 
and the impurity pseudospins. 

It is convenient to describe the staggered component 
of the copper spins by a continuous vector field n (r) 
within the framework of the non- linear ct- model (NLSM). 
Throughout this letter, we use the usual bold font for 
vectors in coordinate space and denote vectors acting in 
spin space by arrows. We use the same orthorhombic co- 
ordinate system in both cases, i.e., Cq. = e^, see Fig. ^ 
Having in mind the above discussion of the relative de- 
grees of freedom, we write the energy of a single layer of 
LSCO as 

1 1 I [Vn (r)]^ + ^ K (r)? + ^ K (r)]' } 
+ V2g / p (r - n) n (r) x k (e^ • V) n (r) , 

(1) 

where n", a = a,b,c, denotes the components of the 
r?- field, subject to the constraint = 1. The first 
line in is the elastic energy which takes into account 
the Dzyaloshinksi-Moriya (DM) and the XY anisotropics, 
with Z? w 2 meV and yT^ « 4 meV. We use ps « 0.18 J 
for the spin stiffness and c « 1.66J for the spin-wave ve- 
locity, with J w 130 meV. These parameters follow from 
neutron scattering data at zero doping jl{j |. The second 
line in represents the interaction of the n-field with 
the impurity pseudospins k — 'I'|(T^'i, see Ref. fl3. The 
impurities are located at positions and p (r) = (r) . 
The coupling constant g k J has been calculated pre- 
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FIG. 1: ( Color online). Characteristic ground state configuration of a particular realization at x = 0.05. The system (a) without 
defects has higher energy than the same system (b) with a vortex-antivortex pair (squares). The impurity pseudospins U are 
oriented along the c axis. Full (open) circles correspond to values —1 (-\-l). Small arrows represent the n-field. The system 
forms domains stretched along the a direction, in which all the pseudospins are aligned in parallel. 



viously within the extended t-J model [l^. Eq. iQJ de- 
scribes the static hmit of the effective field theory. This 
limit is sufficient to study ground state properties be- 
cause ^ results in a long-range interaction [see Eq. (O 
below] . 

The anisotropics in Eq. ^ pin the n-field to a par- 
ticular spatial direction. Because the XY term is larger 
than the DM anisotropy, ^/T^ > D, the staggered field 
n is coplanar and lies in the ab plane. For a single 
Cu02 layer at zero temperature, there is no mecha- 
nism that deflects ft out of this plane. Substituting 
n = (n",n^,n'^) = (sin 0, cos 6*, 0) in Eq. the energy 
of the effective 0(2) NLSM reads 



^ y I [V0 {r f + 2MJ2 hp (r - r,) (e^ • V) (r) | , 

(2) 

where k = ^'fcr'^^'i is an Ising variable taking the values 
±1 and M = V2g/ps « 8. In (0) we have neglected the 
DM term sin^ 9. Despite its smallness, D^/c^ « 10"'', 
this term is important in the Neel phase where it pins the 
n-field to the b axis, i.e., « 0. However, in the spin- 
glass phase, the angle 9{r) quickly varies at a scale much 
shorter than 1]jm — c/D ^ 100 and since (sin^ 9) — 1/2, 
the DM term can be safely neglected. Very importantly, 
Eq. is exact in the sense that it is valid for arbitrary 
6 and not restricted to small angles. 

To integrate out the 6'-field, we start with a single im- 
purity. The variation of Eq. jSJl with respect to 6 yields 
—V'^9 (r) -I- Ml{ef, ■ V)p (r) — 0, which has the solution 



9{y) = I'd (r), where 

^« = ^^{(l + 2-)e-^'"^-l} . (3) 

Because of the linearity of the problem, the solution for 
N impurities is the superposition 

N 

0(r)=^Z,i?(r-r,) . (4) 

1=1 

Substituting this solution into Eq. |5J yields an expres- 
sion of the energy in terms of the Ising pseudospins U 
only 

Ei = ^ hi J {Fi (nrij) + cos (2ay ) F2 (nrij)} , 

(5) 

with Fi (y) = ~y^K2 (2y) and F2 (y) = l/y^-yKs {2y)~ 
y^K2 (2y). Here cky is the angle between the vectors 
and Tij = Ti — Tj, and Kn are modified Bessel functions. 
For r 3> the above expression is equivalent to the 
usual dipole-dipole interaction [cx cos (2aij ) /rfj] . 

Ground state and destruction of Neel order. We per- 
form classical Monte Carlo simulation to find the ground 
state of the Ising pseudospins described by Eq. © us- 
ing the same algorithm as in Re. ^3- Because samples 
of LSCO consist of many Cu02 layers with randomly 
distributed dopants, we average observables over many 
realizations of random impurity distributions. For a par- 
ticular realization, we consider up to iV = 200 Ising pseu- 
dospins h on a square lattice of size L = ^/nJ^. In or- 
der to minimize finite-size effects, we orient the lattice 
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FIG. 2: (a) (n*") as a function of doping x. For k = 0.4, Neel order is destroyed at ~ 0.02, in agreement with experiments. 
For pointlike impurities (k oo), there is no phase transition at least up to x = 0.05. (b) Neutron scattering probability Sq 
for X = 0.04. Full circles correspond to experimental observations taken from Fig. 4 in Ref.\^ with normalized intensities. The 
curve represents our simulation, containing no fitting parameters, (c) Incommensurability 5 (in reciprocal lattice units of the 
tetragonal lattice) as a function of doping. Our calculations (circles) are in good agreement with experimental measurevfients 
(squares) taken from Ref. Q, Q, lid, see Fig. 6 of Ref. Q. Theoretical points for the Neel phase are taken from Ref. 03 



along the orthorhombic coordinate system, apply peri- 
odic boundary conditions along the a axis and extract 
relevant quantities only from the central quarter of the 
system. The h direction is left open, in order not to im- 
pose an artificial constraint on the spiral pitch. Once the 
Ising ground state is found, the n-field is determined ac- 
cording to Eq. Q). Note that 9{y) in can always be 
shifted by a constant Oq. We set 0o = and therefore 
have = and n = Cf, in the undoped system. Fig. ^ia) 
shows a characteristic spin-glass ground state for a par- 
ticular realization at doping x = 0.05. Open (full) circles 
show the Ising pseudospins directed parallel (antiparal- 
lel) to the c axis (orthogonal to the plane). The n-field is 
represented by arrows. The systems forms large domains, 
stretched along the a axis, in which the pseudospins all 
point in the same direction. Fig. ^a) nicely illustrates 
the concept of spiral formation. Despite the random dis- 
tribution of impurities, one easily recognizes horizontal 
lines along which the n-field remains almost constant. In 
contrast to a uniform spiral, these lines are not exactly 
parallel to the a axis, but meander around the impurities. 

The presence of Neel order can be identified as a non- 
zero For a given doping, the average is taken over 
the central quarter of the system and over many realiza- 
tions of random impurity positions. We first performed 
calculations for k — > cx) (pointlike impurities) and found 
that Neel order is not destroyed at least up to a; = 0.05, 
see Fig.|2Ia). In this case, the pseudospins ahgn in an an- 
tiparallel pattern along the h direction, on a scale of the 
order of the separation between impurities. The angle 9 
is thus small and the resulting n-field never completes a 
full rotation. For finite k, the situation is qualitatively 
different because for sufficiently small distances between 
impurities, the combination in brackets {...} in Eq. (O 



is always negative. This favors parallel alignment of the 
pseudospins and hence leads to the formation of parallel 
Ising domains, see Fig. ^a). The domains have a finite 
width because the long-range tail of (|SJ) still favors an- 
tiparallel alignment along the b direction. The width w 
depends on doping and on k, but also on the size of the 
lattice as w (X VL. The domains thus become macro- 
scopic in the thermodynamic limit. The mechanism for 
the domain formation is exactly the same as for ferro- 
magnets |l6j . We will show in the following paragraph 
that topological defects lead to finite, but still very large 
domain sizes. The doping dependence of (n^) for k = 0.4 
is shown in Fig. |2Ia). The phase transition takes place 
at Xc ~ 0.02, in very good agreement with experiments. 
The small deviations from zero at a; > 0.02 in Fig. ^a,) 
are due to finite size effects. It is remarkable that the 
correct critical concentration Xc is obtained for n w 0.4 
that is known from the variable range hopping conduc- 
tivity HI]. 

Topological defects. There are two kinds of topolog- 
ical defects in a 2D 0(3) NLSM, vortices and instan- 
tons [l^. Instantons lead to non-coplanar spin configu- 
rations that, according to our calculations, increase the 
energy because the system prefers coplanar pseudospin 
arrangements. Instantons are therefore not present in the 
ground state. For vortices, the situation is different. Let 
us consider solutions of the Laplace equation V^6' = 
of the form 9^ (r) = Y^jti Qj arg [x - Xj +i{y - Yj)]. 
Here Rj = (Xj , Yj ) and Qj € Z are the positions of the 
impurities and their topological charges (winding num- 
bers), respectively. The energy is minimal for vanishing 
total topological charge with Qj = ±1. In this case, 
we have M/2 pairs of vortices with opposite winding 
numbers. For large average separation between vortices. 
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i? ^ 1, the associated energy calculated with logarithmic 
accuracy is 

Ev = M-Kps In R . (6) 

Because this energy is positive, vortices never appear in 
the ground state of the NLSM without impurities. How- 
ever, in a doped system, one also has to take into account 
the interaction energy Ejv between impurities and vor- 
tices. From Eq. (O, we find 

Eiv = PsM J2 Qj^^tP- [1 - e~'"^" (1 + 2<»,)] , 

i,j=l ^iJ 

(7) 

where — (Rj — ri) are the positions of the vortices 
with respect to the impurities. The interaction energy (jZj) 
is large and negative because of the domains formed in 
the spin-glass phase, in which all pscudospins point in 
the same direction. It therefore favors the creation of 
vortices at zero temperature. The total energy is now 
given by the sum of lO, ® and An illustration 

of a system with one vortex-antivortex pair is shown in 
Fig. mb). Compared to the same realization without 
defects [Fig. QJa)], the energy is around 10% lower. The 
presence of vortices renders the problem of finding the 
ground state quite difhcult. In addition to minimizing 
the energy of the Ising pseudospins, one also has to check 
if vortex-antivortex pairs can lower the energy, and if 
so, optimize their positions. Nevertheless, it is possible 
to estimate the optimal distance between defects in the 
thermodynamic limit. Because of the long-range nature 
of the interaction Q, the energy gain is proportional 
to R. Let us denote the energy of the system without 
vortices by = and use E — Ej + Eiv + Ey for a 
system with M defects. The energy gain due to vortices 
can then be expressed as Ai? = E — E — —-fxMR + 
Mirps InR. By comparing the energies of systems with 4 
and 16 vortices to the same realizations without defects, 
we find that this formula gives the correct R dependence. 
The parameter 7 is weakly doping dependent, at a; = 
0.02, we obtain 7 « 2.0, leading to an optimal distance 
Ropt ~ 130 in the thermodynamic limit. At x — 0.05, 
we find 7 « 1.1 from which we deduce i?opt ~ 90. We 
thus conclude that vortex-antivortex pairs break up the 
domains and lead to a parquet-like arrangement of Ising 
pseudospins with a natural domain size R ~ 100. 

Static spin structure factor and neutron scattering. Let 
us finally calculate the structure factor 

as well as the sum S'q = SapS^^, which is related to the 
neutron scattering cross section. Fig^^E^b) shows exper- 
imental data taken by Fujita et al. [j| on a sample with 
doping X = 0.04 together with our theoretical results for 



S'q. The agreement between theory and experiment is 
quite remarkable, especially given the fact that it is not 
restricted to a particular sample, but can be observed 
over a broad range of doping. This latter finding is illus- 
trated in Fig.|2Ic), which shows the incommensurability 
(5, defined as half the distance between the peaks, as a 
function of doping. Expressed in reciprocal lattice units 
of the tetragonal structure, the incommensurability is in 
good approximation proportional to doping, (5 ~ x. Be- 
cause of their low density, topological defects do not have 
a major influence on S'q, and only lead to some broaden- 
ing of the peaks. 

To conclude, we have developed a low energy effec- 
tive field theory that describes the magnetic structure of 
La2_2:Sr2,Cu04 in the spin-glass phase (x < 0.055). We 
have shown that the staggered component of the cop- 
per spins is confined to the ab plane, due to the XY 
anisotropy. The static spin structure factors obtained 
in our approach are in excellent agreement with neu- 
tron scattering data, over a whole range of doping. We 
have analyzed the nature of the quantum phase transi- 
tion from the Neel to the spin-glass phase and have also 
shown that topological defects (spin vortices) play a sig- 
nificant role in the ground state of the spin-glass phase. 
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